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ABSTRACT 

At redshifts z < 1 a large fraction of the baryons is thought to reside in diffuse gas that has 
been heated to temperatures T ~ 10^ - 10*" K by gravitational accretion shocks and galactic 
winds. Absorption by the 770.41 A, 780.32 A doublet of Ne vm in the spectra of quasars 
provides us with a unique tool to study this elusive gas phase. We construct an analytic model 
for the Ne vm absorption as a function of the gas density, temperature, and the abundance 
of heavy elements, allowing for the metals to be in-homogeneously distributed within the 
gaseous structures responsible for the absorption, and assuming that these structures are self- 
gravitating, that the gas is exposed to the UV/X-ray background, and that it is in ionisation 
equiUbrium. Our analytic model agrees with the predictions of a cosmological, hydrodynamic 
simulation from the Overwhelmingly Large Simulations (OWLS) project provided that the 
degree of metal mixing, parametrized by the average Une-of-sight metal-filling fraction within 
the gas clouds harbouring Neviii, is low: ci ~ 10 Most of the Nevm at low redshift is 
produced by collisional ionisation in gas with densities 10"'' < (riu/cm^^) < 10"^, temperatures 
10^ < (T/K) < 10^, and column densities Nnc VIII < 10 cm . Absorbers as strong as those 
recently detected by the Cosmic Origins Spectrograph on HST, NNevm ^ lO'"* cm"^ (or rest- 
frame equivalent widths Wr > 30 mA), are only produced in collisionally ionised gas with 
temperatures T 5 x 10^ K and tin Z 10""* cm"-'. The neon abundance of the Ne vm absorbers 
in our simulation is distributed over the range - 1 .5 < [Ne/H] < 0, with median [Ne/H] = -0.7. 
Metal-rich ([Ne/H] « 0), strong (NNe™ ~ 10^* cm"^) absorbers have typical sizes ~ 10 kpc 
and likely originate in the immediate vicinity of galaxies. Despite their relatively high total 
hydrogen column densities (A'h ~ 10^*^ cm"^), these strong Nevm cloudlets account for less 
than 1 per cent of the cosmic baryon budget. The baryon content of the gas harbouring the 
Ne VIII absorbers (and which has a similar density and temperature as the Ne viii cloudlets) is, 
however, a factor higher. We conclude that strong Ne vm absorbers are robust probes of 
shock-heated diffuse gas. Spectra with signal-to-noise ratio ~ 10^ would allow the detection of 
the weaker systems that trace a substantial fraction of the baryons at low redshift. 

Key words: cosmology: theory — methods: numerical — analytical — quasars: absorption 
lines — galaxies: formation — intergalactic medium 



1 INTRODUCTION 

The analysis of the physical conditions and the distribution and 
abundances of heavy elements in interstellar- and intergalactic gas 
using absorption spectroscopy is a well-established technique which 
has yielded a wealth of information about the status and the evolution 
of a significant fraction of the baryonic matter in the Universe. 
The Cosmic Origins Spectrograph (COS; Green & Morse 1998; 
Green et al. 2012), an ultra-violet (UV) spectrograph installed on the 
Hubble Space Telescope (HST), has opened a new era of absorption 
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spectroscopy, as a result of its increased sensitivity (from factors of 
a few to an order of magnitude) with respect to previous instruments 
such as the Far Ultraviolet Spectroscopic Explorer (FUSE) satellite 
or the Space Telescope Imaging Spectrograph (STIS) installed on 
HST. COS is designed to perform high-sensitivity, medium- and 
low-resolution (FWHM a; 1 7 km s" ' ) spectroscopy of astronoinical 
objects in the 1 150 A - 3200 A wavelength range. 

One of the major science goals of COS is to clarify the long 
standing mystery of the whereabouts of a significant fraction (30 
to 50 per cent) of the baryons in the low-redshift Universe (Persic 
& Salucci 1992; Fukugita 2004; Shull et al. 201 1), that were syn- 
thesised in the Big Bang and that are detected mainly through Hi 
absorption at redshift z ^ 3 (Rauch et al. 1997; Weinberg et al. 1997; 
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Schaye 2001). A robust prediction of cosmological simulations is 
that most of these "missing baryons" reside in a gas phase with 
densities nn ~ 10"* - 10"'cm~^ that, as a result of the formation 
of structure in the Universe, has been shock-heated to temperatures 
r > 10' K (e.g. Cen & Ostriker 1999; Dave et al. 2001; Bertone et al. 
2008). Along with gravitational shock-heating, feedback processes 
such as supernova (SN) explosions and outflows driven by active 
galactic nuclei (AGN) have been shown to significantly contribute 
to the mass in this difl\ise, warm-hot gas phase at low redshift (e.g. 
Tepper-Garcfa et al. 2012). 

So far, attempts to trace shock-heated diffuse gas at low red- 
shift, and to constrain its baryon content, have mainly relied on 
the detection of absorption by five-times ionised oxygen (O vi) in 
the UV spectra of distant point-like sources (e.g. quasars, QSOs; 
Tripp et al. 2000; Richter et al. 2004; Danforth et al. 2006; Danforth 
& Shull 2008; Thom & Chen 2008b; Tripp et al. 2008; Danforth 
et al. 2010; Prochaska et al. 201 1). However, because O vi may be 
produced both by photo-ionisation in gas at T ~ 10"* K and by col- 
lisional ionisation in gas at 7 a; 3 x 10' K, and because the metals 
may well be poorly mixed on small scales (Schaye et al. 2007), 
the interpretation of intervening O vi absorbers and their status as 
tracers of shock-heated gas is still controversial, as shown by both 
observational (Thom & Chen 2008a; Tripp et al. 2008; Danforth 
&. Shull 2008) and theoretical studies (e.g. Oppenheimer & Dave 
2009; Tepper-Garci'a et al. 201 1 ; Smith et al. 201 1 ; Cen 2012). Be- 
sides Ovi, broad {b^i > 40km s"') and shallow (to(Hi) <, 0.1) 
Lyor absorption features, the so-called Broad Lya Absorbers (BLAs; 
Richter et al. 2006), have been considered as a means to detect 
diffuse warm-hot gas at low redshift (Sembach et al. 2004; Richter 
et al. 2004; Williger et al. 2006; Lehner et al. 2007; Danforth et al. 
2010; Tepper-Garcfa et al. 2012). However, BLAs are difficult to 
identify in the UV spectra of QSOs because of the limited signal- 
to-noise (S/N) and spectral resolution of data obtained with current 
space-based UV spectrographs. 

An alternative to O vi and BLAs as tracers of shock-heated 
gas might be offered by seven-times ionised neon (Ne vm). Neon is 
the fifth most abundant element in the Universe,' and this particu- 
lar ionisation state has two resonant transitions with rest-frame 
wavelengths in the extreme ultra-violet (EUV) at wavelengths 
770.41 A and 780.32 A (e.g. Verner et al. 1994), which are red- 
shifted into the COS spectral range at z > 0.5. The Ne vm ionisation 
fraction in coUisional ionisation equilibrium (CIE) is highest at 
temperatures T « 5 x 10' K (e.g. Sutherland & Dopita 1993). 

To date, though, only a handful of intervening (as opposed to 
'proximate to the QSO') Ne vm absorbers have been identified in UV 
absorption spectra at redshifts z < 1 (Savage et al. 2005; Narayanan 
et al. 2009, 201 1 ; Tripp et al. 201 1; Meiring et al. 2012). Interest- 
ingly, all these absorbers have similar strengths, with column densit- 
ies NNevm ~ 10'*cm~^, and they all show a simple kinematic struc- 
ture (at the given sensitivity and spectral resolution), thus suggesting 
that the observed Ne vm absorption is produced in gas under rather 
uniform physical conditions (but see Tripp et al. 201 1). In particular, 
the narrow range of measured temperatures (logiQ(r/ K) « 5.4-5.7) 
indicates, first, that Ne vm absorbers are very sensitive probes of the 
thermal state of the absorbing gas; and second, that Ne vm absorbers 
are reliable tracers of gas at T « 5 x 10' K (Savage et al. 2005, 
but see Oppenheimer & Schaye 2013a who argue that Nevm ab- 
sorbers may trace fossil AGN proximity zones). Furthermore, their 
relatively high hydrogen column densities (A^h ~ - lO^^cm"^) 
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suggest that they contain a substantial amount of baryons (Savage 
et al. 2005; Narayanan et al. 2009). In view of these findings, Ne vm 
absorbers indeed seem a promising tool in the search for warm-hot, 
diffuse gas and the 'missing baryons' in the low redshift Universe. 

A clear link between the intervening Ne vm absorbers and this 
peculiar gas phase has, however, not as yet been observationally es- 
tablished due to the current small number of detections, and support 
from the theoretical side is hence called for. This is the motivation 
for the present paper. 

Here, we investigate in detail the physical conditions of gas 
producing Ne vm absorption at low redshift following two rather 
independent and complementary approaches. First, we make use 
of a cosmological simulation that samples a representative volume 
of the universe at z = 0.5 and which includes many of the physical 
processes believed to be important for the production and distribu- 
tion of metals in the Universe. From this simulation we extract the 
distribution of densities and temperatures, the baryon content and 
the level of enrichment of a statistically significant sample of Nevm 
absorbers. Additionally, we use this sample to make predictions for 
two canonical absorption-line statistical quantities: the distribution 
of Ne vm column densities and the number of Ne vm absorbers of 
a given strength per unit redshift. Second, we develop an analytic 
model to predict standard observables (i.e. Nevm column density, 
Nevni central optical depth tq, rest-frame equivalent width W^) 
of Ne vm absorption features as a function of the density and the 
temperature of the absorbing gas. The fundamental assumptions 
underlying our analytic model are that the gas phase producing 
Nevm absorption: 1) is exposed only to the extragalactic UV back- 
ground; 2) is in ionisation equilibrium; 3) is embedded within a 
larger gaseous structure in (quasi) local hydrostatic equilibrium; 4) 
has an inhomogeneous metal distribution, with the metals filling only 
a fraction of the volume projected along the sightline across the 
structure. We present and discuss these two approaches in Sections 
2.1 and 2.2, respectively. 

The combined results from our simulation and our analytic 
model -which we discuss in Section 3- allow us to determine the 
gas phases in the low redshift Universe that are expected to contain 
significant amounts of Ne vm (in term of mass), and to identify these 
gas phases through their Ne vm absorption. Our ultimate goal is to 
aid the interpretation of current observational data (discussed in 
Section 4), and to provide a guide for future observations aiming to 
detect Ne vm absorbing gas using absorption spectroscopy. 

We refer the impatient reader to Section 5, where we summarise 
the central results of our study. A detailed presentation of the status 
of current Ne vm detections at low redshift and other results that 
might be of interest only to the more speciaUsed reader are left for 
the appendix. 



2 LOW-REDSHIFT Ne vm ABSORBING GAS: 
SIMULATION AND ANALYTIC MODEL 

Seven-times ionised neon (Ne vm) can be produced by the ionisa- 
tion of Nevn (ionisation potential Ep = 207.271 ± 0.013 eV; e.g. 
Kramida & Buchet-Poulizac 2006) via the absorption of photons 
with energies Ey > Ep and by collisions with electrons and other 
atoms. The Nevm ion fraction /Nevm = (nNevm/wNe) produced by 
both photo-ionisation and coUisional ionisation in gas as a function 
of the gas density and temperature is presented in Figure 1, which 
we obtain from the tables provided by Wiersma et al. (2009a). These 
tables, which were also used to calculate the ionisation balance in 
our simulation (see Section 2.1), were computed with the photoion- 



© — RAS, MNRAS 000, 1-17 



Ne vm absorbers 3 




Figure 1. Nevm ionisation fraction, ^evm = («Neviii/nNc), of optically 
tliin gas in ionisation equilibrium exposed to the X-Ray/UV background 
radiation from galaxies and quasars at z = 0.5. The numbers in the panel 
indicate the minimum /Nevm (in logarithmic scale) in gas at the densities 
and temperatures within the corresponding contour Note that relatively high 
ionisation fractions /nc viii ~ 0. 1 occur in both low-density, low-temperature 
(photo-ionised) gas with mh < 10"^ cm"-' and T < 10^ K, and in denser, 
warm-hot (colUsionaUy ionised) gas with mh > 10 cm" -3 andr K SxlO^K. 

isation package cloudy (version 07.02.00 of the code last described 
by Ferland et al. 1998) assuming the gas is exposed to the z = 0.5 
Haardt & Madau (2001) model for the X-Ray/UV background ra- 
diation from quasars and galaxies (assuming a 10 per cent escape 
fraction for H-ionising photons). In addition, the gas is assumed to 
be dust-free, optically thin and in ionisation equilibrium. The most 
striking feature visible in Figure 1 is the high Ne vm ion fractions, 
/Nevm ~ 0.1, both in photo-ionised gas with nn < 10"' cm"^ and 
T < 10' K, and in collisionally ionised gas with > 10"' cm"^ 
and r » 5 X 10' K. This is remarkable because it indicates that the 
presence of measurable Ne vm absorption in spectra may by itself 
not be a sufficient condition for the detection of gas at T > 10' K. 

It is therefore important to explore in detail the connection 
between the physical conditions of the gas harbouring Ne vin and 
the characteristics of the Ne viii absorption it produces, which we 
do next using numerical simulations and analytic methods. 

2.1 Simulation 

The simulation used in this work is one of the several cosmological 
simulations that together comprise the Overwhelmingly Large Sim- 
ulations (OWLS) project, described in detail in Schaye et al. (2010). 
These simulations were performed with a significantly extended ver- 
sion of the Af-Body, Tree-PM, Smoothed Particle Hydrodynamics 
(SPH) code GADGET in - which is a modified version of gadget n (last 
described in Springel 2005) -, a Lagrangian code used to calculate 
gravitational and hydrodynamic forces on a system of particles. 

Here we focus on a particular model from the OWLS suite 
referred to as AGN. This model adopts a flat ACDM cosmology 
characterised by the set of parameters [Q.^, fib, f^A, erg, Ws, h] = 
(0.238, 0.0418, 0.762, 0.74, 0.95, 0.73) as derived from the Wilkin- 
son Microwave Anisotropy Probe (WMAP) 3-year data^ (Spergel 

^ These parameter values are largely consistent with the WMAP 7-year 



et al. 2007). This model includes star formation following Schaye & 
Dalla Vecchia (2008), metal production and timed release of mass 
and heavy elements by intermediate mass stars, i.e., asymptotic 
giant-branch (AGB) stars and supernovae of Type la (SNIa), and 
by massive stars (core-collapse supernovae, SNIIe) as described by 
Wiersma et al. (2009b). It further incorporates kinetic feedback by 
SNIIe based on the method of Dalla Vecchia & Schaye (2008), as 
well as thermal feedback by SNIa (Wiersma et al. 2009b). Radiative 
cooling by hydrogen, helium and heavy elements is included ele- 
ment by element following the method of Wiersma et al. (2009a). 
The ionisation balance for each SPH particle is computed as a func- 
tion of redshift, density, and temperature using pre-computed tables 
obtained with the photoionisation package cloudy (version 07.02.00 
of the code last described by Ferland et al. 1998), assuming the gas 
is exposed to the Haardt & Madau (2001) model for the X-Ray/UV 
background radiation from galaxies and quasars. Additionally, our 
simulation includes feedback by active galactic nuclei (AGN) based 
on the model of black hole growth developed by Booth & Schaye 
(2009), which is a modified version of the model by Springel et al. 
(2005). Our simulation was run in a cubic box of 100/z"' comoving 
mega-parsec (Mpc) on a side, containing 512^ dark matter (DM) 
particles and equally many baryonic particles. The initial particle 
mass is 4.1 x 10*/!"' Mq (DM) and 8.7 x lO'^/i"' Mq (baryonic). The 
gravitational softening is set to 8 /?"' comoving kUo-parsec (kpc) 
and is fixed at 2/;"' proper kpc below z = 3. 

It is worth noting that this model has been shown to self- 
consistently reproduce several fundamental measurements, namely: 
the observed mass density in black holes at z = 0; the black hole 
scaling relations (Booth & Schaye 2009) and their evolution (Booth 
& Schaye 201 1); the observed optical and X-ray properties, stellar- 
mass fractions, star-formation rates (SFRs), stellar-age distribu- 
tions and the thermodynamic profiles of groups of galaxies (Mc- 
Carthy et al. 2010); and the standard Hi observables (distribution 
of Hi column densities A'hi, distribution of Doppler parameters bm, 
bm - Nm correlation) and the BLA line number density at low red- 
shift (Tepper-Garcfa et al. 2012). For a more detailed description 
of this (and other) models that are part of the OWLS project, see 
Schaye et al (2010). 

2.1.1 Physical conditions ofNe vm absorbing gas 

We investigate the physical state of the gas containing Ne vm in 
our simulation at z = 0.5 in two different ways. Our first method is 
to compute the distribution of Ne viii mass in gas as a function of 
density and temperature using the density, temperature, and Ne vm 
mass of each SPH particle. For the second method we draw a number 
(5000) of random sightlines through our simulation box and generate 
and fit a corresponding set of high-S/N, COS-like spectra, using the 
procedure described in Appendix B. In this way we obtain a sample 
of roughly 8670 Ne vm absorption components and we assign each 
a characteristic (i.e. Nevm optical depth weighted ) density (hh) 
and a characteristic temperature (T). To assign a Ne vm optical 
depth weighted density (temperature) to an absorber, we follow the 
procedure described in Tepper-Garcfa et al. (2011, their Section 
5.1). Briefly, we use our synthetic spectra to compute the density 
(temperature) weighted by the Ne vm optical depth in redshift space 
along the sightline as in Schaye et al. (1999). We then identify 

results (Jarosik et al. 2011), the largest difference being the value of trg, 
which is 2 cr lower in the WMAP 3-year data than allowed by the WMAP 
7-year data. 
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Figure 2. Top-left: Distribution of Ne viii mass in gas in our simulation at z = 0.5. The coloured regions show the fraction of Ne vin mass in gas contained within 
the corresponding region of im - T space, as indicated by the colour bar in the bottom right margin. The soUd contours show the distribution (by number) of gas 
densities and temperatures of the Ne vin absorbers detected in a set of synthetic spectra at z = 0.5; the outer-most contour contains 90 per cent of the absorbers. 
The blue dashed curves indicate the minimum column density logjg(NNeviii/cm"^) = 11 allowed in our fitting procedure. Note that virtually all the Ne viii mass 
in gas and all the detected Ne viii absorbers are contained within the region defined by this limit. Top-right: Scatter plot of Ne viii absorbers colour-coded by their 
Ne VIII column density. The solid contours indicate the Ne viii column density (in cm~^) predicted by our analytic model, equation (7), adopting [Ne/H] = -0.7 
(relative to logm («Nc/'Ih)o = -3.91), fg = 0. 168 and a filling fraction logjQ ct = -0.5. Bottom-left: Scatter plot of Ne viii absorbers colour-coded by their 
central optical depth. The contours show the Ne viii central optical depth predicted by our analytic model, equation (10; the values for [Ne/H], fg and are 
identical to the values adopted before). Bottom-right: Scatter plot of Ne viii absorbers colour-coded by their rest-frame equivalent width (in mA). The contours 
indicate the Ne viii rest-frame equivalent width in mA computed from the Ne viii column density predicted by our analytic model, equation (11). Note that the 
Ne VIII gas mass distribution in the top-left panel, and the contours in each panel are identical to the coiTesponding contours shown in Figure 8. The contours 
agree with the coloured points, indicating that our analytic model reproduces the simulation results very well. 



the pixel corresponding to the absorption line's centre, smooth (i.e. 
average) over ±2 pixels around this pixel, and assign this last optical 
depth weighted average density (temperature) to the absorber. Note 
that our synthetic spectra have a constant pixel size Sv = 3 km s"' , 
which corresponds to a physical linear path SI « 33 kpc at z = 
0.5, assuming that the velocity spread is due solely to the Hubble 
expansion (note that our spectra are computed taking into account 
also peculiar velocities and thermal broadening). We thus smooth 
over a range Av = 12 km s"' or 132 kpc around the line centre. 

The distributions of densities (na) and temperatures (T) of the 
Ne VIII gas in our simulation at z = 0.5 obtained from these two 
different approaches are compared in the top-left panel of Figure 2 
(note that the other panels will be discussed in Section 2.2). The dis- 
tribution of Ne VIII mass as a function of and T obtained from the 
density, the temperature, and the Ne viii mass of each SPH particle 
is shown by the coloured areas. Each of these is colour coded by the 



fraction of Ne viii mass in gas contained within the corresponding 
region of na- T space, according to the colour bar displayed in the 
bottom right margin. So, for example, the innermost region (orange) 
contains 25 per cent of the Ne viii mass in the simulation; the cyan 
and the orange area together enclose 50 per cent of the Ne vm mass, 
and so on. The black contours show the distribution of Ne vm optical 
depth weighted densities and temperatures of the Ne viii absorbers 
detected in our simulated spectra. These contours enclose, starting 
from the innermost, 20, 40, 60, 80, and 90 per cent of the total 
number of absorbers. Note that the physical conditions of the Ne vm 
absorbing gas inferred from the spectra (black contours) match the 
corresponding intrinsic properties of the Ne viii harbouring gas in 
the simulation (coloured regions) very well. This confirms that our 
choice of assigning physical quantities to the absorbers identified 
in synthetic spectra by using optical depth weighted quantities is a 
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Figure 3. Comparison of the density (left), temperature (centre), and tliermal pressure (right) of the Ne vin absorbing gas and that of the gas responsible for the 
H I absorption at the same location in redshift space. 



reliable method to unequivocally link the absorption of the gas to its 
physical state. 

We find that most of the Ne vm mass in our simulation resides 
in gas with densities with densities 10"* < (nn/cm"^) < 10""* (or 
roughly 1-150 times the cosmic mean at z = 0.5) and temperatures 
10' % (T/K) < 10*. Hence, even though the Nevm ion fraction 
can be high both in photo- and coUisionally ionised gas (see Figure 
1), our simulation predicts that most of the Nevm is produced by 
coUisional ionisation in low-density gas in ionisation equilibrium. 
This result disagrees in part with Oppenheimer et al. (2012), who 
find that most of the Ne vm mass in their simulation (albeit at z = 1) 
is also in gas with densities Mh ~ 10"' cm"', but at significantly 
lower temperatures, only slightly above T ~ i(f K. We note that 
we do find in our simulation Ne vni absorbers with T < 10' K (see 
scatter plots in Figure 2); however, these are generally weak and 
their relative number is low. 

A similar discrepancy concerning the typical temperature of 
O VI absorbing gas between Oppenheimer et al.'s simulations and 
ours was reported in Tepper-Garcia et al. (201 1). Oppenheimer et al. 
attribute the low fraction of highly-ionised metal absorbers with 
r ^ 10' K in their simulation to the fact that their simulation does 
not resolve conductive hitcrfaccs within hales. Since our simulation 
does not resolve such interfaces either, the difference in the predicted 
temperature distribution of Ne vm (or O vi) absorbing gas most 
probably stems from the difl^erence in the implementations of the 
sub-grid physics; in particular, the prescription to model winds on 
galactic scales may be crucial since these determine the distribution 
of mass and heavy elements among the different gas phases. Also, 
the assumption of ionisation equilibrium, which is also used by 
Oppenheimer et al. (2012), may be important (e.g. Oppenheimer 
& Schaye 2013b,a) and future studies relaxing this assumption 
may obtain different results. We therefore advise that the physical 
conditions, in particular the thermal state, of highly ionised gas 
predicted by simulations should be taken with caution. 

As will be explained in more detail in Section 2.2, two of the 
implicit assumptions underlying our analytic model are that the 
density («h) and the temperature {T) of the Ne vni absorbing clouds 
are comparable to the local density and the local temperature of the 
surrounding gas. This is reflected in the fact that the Ne vm column 
density (and all the other derived spectral properties) as given by 
equation (7) depends on n^Lj, where Lj itself depends on wh and 
T (see equation 4). Moreover, since rinLj oc {nnTy^ = {P/ky^ 
(where P is the thermal pressure and k is Boltzmarm's constant), 
these assumptions imply that the metal absorbers are not far from 
thermal pressure equiUbrium with their surrounding (i.e. confining) 
medium. 

In order to investigate if these conditions are satisfied by the 



Ne vm absorbing gas in our simulation, we compare the character- 
istic density (wh). temperature (T), and thermal pressure {P/k = 
nuT) of each Nevm absorber in our sample with the correspond- 
ing quantity of the associated (i.e. aligned in velocity space) Hi 
absorber. To this end, we compute the Ne vm- and the H i-optical 
depth weighted density and temperature at the centre of each Nevm 
absorbing component identified in our simulated spectra, following 
the previously described procedure. Note that when comparing Hi 
to Ne vm weighted quantities, it is implicitly assumed that relative 
peculiar velocities between the Nevm and Hi absorbing phases 
are small, so that their alignment in velocity space implies their 
coincidence in physical space. The comparison of the densities, 
temperatures, and the thermal pressures thus obtained is shown in 
Figure 3. In each panel, the x-axis indicates the physical quantity 
(e.g. density) weighted by the Ne vm optical depth and the y-axis in- 
dicates the corresponding quantity weighted by the Hi optical depth. 
The contours lines show the distribution (by number) of absorbers; 
these contours enclose, starting from the innermost, 20, 40, 60, 80, 
and 90 per cent of the total number of absorbers. The dotted line 
along the diagonal indicates in each case a perfect match between 
the two physical quantities considered. 

The left panel of Figure 3 shows that the density of the Ne vm 
absorbing gas is similar to the density of the surrounding gas for 
most of the Ne vin absorbers, being only slightly higher on average 
in the Hi absorbing phase with nn ^ 10"' cm"'. In the middle 
panel we see that, while the temperature of both these phases is 
comparable for most of the absorbers, for a non-negligible fraction 
of the absorbers it is a factor ~ 10 higher in the Nevm absorbing 
phase than in the H i absorbing phase. The right panel demonstrates 
that the thermal pressure (P/k) is very similar in both gas phases; 
in other words, the phases are, on average, in pressure equilibrium 
with each other. These results suggest that a typical Ne vm absorber 
corresponds to a metal-emiched pocket of warm-hot gas confined 
by a slightly denser and cooler gaseous structure. 



2.1.2 Absorption line statistics 

Using the sample of Ne vni absorption components obtained from 
our set of synthetic spectra (see Section 2.1 and Appendix B), 
we compute the Nevm line number density and the distribution 
of Ne VIII column densities at z = 0.5 . The results are shown 
in Figures 5 and 6, respectively, and discussed below. An ex- 
ample of a strong Ne vm absorption feature produced in gas with 
logiQ(nH/cm"') s! -3.5 and logiQ(7'/K) » 6.1 along a random 
sightline in our simulation at z = 0.5 and its decomposition into 
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Figure 5. Number of Ne viii absorption components per unit redshift above 
a given rest-frame equivalent width (dN/dz{W > W, ); black curve) and 
the corresponding Poisson uncertainty (shaded area) obtained from a set of 
simulated COS-like spectra with S/N=1000 along 5000 random sightlines 
drawn from our simulation atz = 0.5. The top ,v-axis indicates the approx- 
imate Ne vni column density corresponding to a given rest-frame equivalent 
width. The blue dashed curve shows the probability distribution of Ne vni 
equivalent widths; its amplitude is indicated by the y-axis on the right. The 
vertical dotted line indicates the median value Wr = 1 .23 mA. Note the small 
fraction of absorbers with > 10 mA. 



Figure 4. Nevin absorption feature with total log[o(NNcviii/cm" ) « 14.4 
along a random sightline drawn from our simulation at z = 0.5 in a synthetic 
COS-like spectrum adopting different S/N values (black histogram) and 
the corresponding fits (red). The vertical dashes mark the velocity centre 
of each fitted component; their corresponding logarithmic column density, 
rest-frame equivalent width (in mA), and Doppler parameter (in kms"') are 
listed in the bottom-left of each panel. While the equivalent width (W^) and 
column density (NNcvm) integrated across the absorber are roughly the same 
for dift'erent S/N values (Wr » 105 mA and log,(|(NNcviii/cm"^) k 14.4, 
respectively), the values of Wr and Nncviii of the individual components vary 
substantially from panel to panel since a different number of components is 
required in each case to produce a good fit. Note that all fits have a reduced 
;f ~ 1 (see Appendix B for details on the fitting procedure). This peculiar 
absoiption feature is produced in gas with log[Q(nH/cm"^) sb -3.5 and 
logio(r/K)^6.1. 

individual Gaussian components is shown in Figure 4 (see the figure 
caption for more details). 

Figure 5 displays the number of Ne vin absorption components 
per unit redshift with rest-frame equivalent width above a given 
threshold, dN/dz(W > W,-), predicted by our simulation at z = 0.5 
(black) and the corresponding Poisson uncertainty (gray shaded 
area). The blue dashed curve shows the probability distribution 
(indicated by the right y-axis) of rest-frame equivalent widths, which 
has a median value Wr = 1.23 mA (indicated by the vertical dashed 
line). Thus, roughly half of the Nevni absorbers in our high-S/N 
(i.e. S/N=1000) simulated spectra have IV,. < 1 mA, while the 
other half have a rest-frame equivalent width in the range 1 mA < 
Wr i 10 mA, with a very small fraction having IV^ > 10 mA, 
dN/dz(W > 10 mA) ~ 1, and an even smaller fraction with Wr > 
30 mA, dN/dziW > 30 mA) ~ 10^'. 



The value of the cumulative line frequency for > 30 mA 
we find is consistent with the predictions from other simulations 
(Oppenheimer et al. 2012), but significantly lower than the values 
obtained from current observations. Narayanan et al. (2009) give 
dN/dz(W > 30 mA) = 2.1 (no error quoted), which they estimate 
indirectly from the Ovi line-number density at comparable red- 
shifts, assuming every seventh O vi absorber has associated Ne viii 
absorption. Based on three detections along a single sightline, Meir- 
ing et al. (2012) obtain dN/dz(W > 30 mA) = 7+^, which is only 
consistent with the value given by Narayanan et al. (2009) at the 
lower 2cr level. The poor agreement between these two measure- 
ments stems most probably from the fact that they rely on very few 
detections along small total redshift paths, and they are therefore 
highly uncertain. In particular, it is possible that the measured line 
frequencies have been overestimated. Increasing the sensitivity of 
the data would reveal weaker Ne viii absorbers along many more 
sightlines, thus increasing the number of detections but also the 
surveyed redshift path. It is, however, not clear whether the meas- 
ured line frequency would converge to a value which is more in line 
with the result from our simulation. Alternatively, we consider the 
possibility that the low rate of Ne viii absorbers with Wr ^ 30 mA 
(or log(Q(NNcviii/cm"-) > 14) in our simulation might be an artefact 
related to the lack of high equivalent width O vi absorbers in our 
(Tepper-Garcfa et al. 201 1) and other (e.g. Oppenheimer & Dave 
2009) simulations that may result from turbulent mechanisms on 
sub-resolution scales which are not properly captured by these sim- 
ulations. It is also possible that the observed strong Nevin (and 
O vi) absorbers trace fossil AGN proximity zones (Oppenheimer & 
Schaye 2013a). Ultimately, however, a significantly larger volume 
of data is required to get a reliable estimate of dN jdziW > W, ), and 
to test our predicted value of dN/dz(W > 30 mA) ~ 10"' . 

Figure 6 shows the distribution of Neviii column densit- 
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Figure 6. Distribution of Ne viii column densities per unit Nnc viii per unit 
absorption distance, /(NNeviii) = (i'N ldx(i^ticw\\\ of the components iden- 
tified in a set of simulated high-S/N. COS-like spectra along 5000 random 
sightlines drawn from our simulation at z = 0.5 (red symbols). The black 
symbols show the distribution of column densities integrated (i.e. adding 
up all the fitted components) along each spectrum; the blue curve show the 
distribution of intrinsic column densities integrated along each correspond- 
ing physical sightline; the error bars along the y-axis and the shaded area 
indicate the corresponding Poisson scatter. The dashed line shows a fit to the 
red symbols in the range logjQ(NNevni/cm"^) £ [12.5, 14.5] in the form of 
a power-law with slope yS = -2.9 ±0.1. Note that the line has been shifted 
upwards for display purposes. 



the black and the blue distributions would indicate that the total- 
ity of the Ne vm (mass) is detected in absorption. The recovery is 
excellent for NNevm > lO'^cm"^. However, while the distribution 
of integrated, intrinsic column densities (blue / shaded) flattens off 
at column densities Nncviii ~ lO'^cm"^, the distribution of column 
densities measured from the spectra (black / red) drops sharply' at 
Nncviii :£ lO'^cm"'. Thus, Ne VIII absorbers with Nfje VIII :£ lO'^cm"^ 
should exist but their detection in absorption spectra requires S/N > 
1000. From the baryon content in the form of Ne viii in the simula- 
tion, ^^"y,,, = 2.70 X 10"**, and the baryon content in the form of 
Ne VIII recovered from the spectra, ^^"^^1 = 2.25x 10"**, we estimate 
that the fraction of Ne VIII mass in absorbers with Nfje VIII < lO'^cm"^ 
is roughly 20 per cent. 



2.2 Analytic model 

Now that we have discussed the physical conditions of the gas 
containing Ne viii in our simulation, we proceed to investigate the 
connection between the physical state and the absorption signatures 
of this gas phase using an independent approach. As we have done 
in Tepper-Garcfa et al. (2012) for H i, our goal here is to develop a 
model that allows us to predict the Ne viii line observables (column 
density, equivalent width, central optical depth), and the spectral 
sensitivity in terms of S/N required to detect Neviii absorbing 
gas as a function of the gas density («h) and the gas temperature 
(r). Conversely, such a model should allow one to constrain the 
physical state of the absorbing gas from its observed absorption 
characteristics. 



ies per unit NNevm per unit absorption distance,^^ /(Nncviii) = 
(5^A/'/o(y5NNeviii, in our simulation at z = 0.5. The red symbols 
show the distribution of Ne vm column densities (CDDF) of the 
components identified in a set of high-S/N, COS-like spectra along 
5000 random sightlines. The dashed line corresponds to the fit in the 
form of a power-law with slope /3 = -2.9 ± 0. 1 to the red symbols in 
the range log[(,(Ne vm/cm"^) 6 [12.5, 14.5], and it has been shifted 
upwards for display purposes. We note that the slope /? = -2.9 ±0.1 
of the Ne vm CDDF at z = 0.5 is significantly steeper"* than the 
slope p = -2. 1 ± 0. 1 of the O vi CDDF obtained from observations 
(Tilton et al. 2012) and found in our simulations at a somewhat 
lower redshift (z = 0.25; Tepper-Garcfa et al. 2011). 

The black symbols in Figure 6 show the distribution of column 
densities integrated (i.e. adding up the column densities of all identi- 
fied components) along each spectrum of size 8340 km s"' atz = 0.5. 
The error bars along the y-axis, indicate the corresponding Poisson 
uncertainty. The blue curve and the shaded area show the distri- 
bution of intrinsic Ne vni column densities integrated along each 
corresponding physical sightline of length 100/j"' Mpc (comoving) 
and its Poisson uncertainty, respectively. A perfect match between 

^ The absorption distance in a flat Universe dominated by matter and Dark 
Energy is given by 

* = (1 + z)\ci„,(\ + zf + nA]-'/2rfz . 

^ The slope found by Tilton et al. (2012) and Tepper-Garcfa et al. (201 1) is 
obtained from the analysis of observed and simulated spectra, respectively, 
with moderate S/N; 10-30. Using the set of high-S/N, COS-like synthetic 
spectra of this study, we find /J = -2.0 ±0.1 for the Ovi CDDF, which 
is consistent with the slope obtained in our previous work from simulated 
spectra with S/N = 10 - 30. 



The Nevm central optical depth of an absorption line with 
column density NNevm and Doppler parameter ^Neviii is given by 

/XT ^ /- 1 ^ /Nnc villi 

To(Ne vm) = ^jn r, c /o 

\ ONeVIII / 

^. n ^ne. i ^Neviii \{ ^Neviii '\ 

where r^ = 2.82 x 10"" cm is the classical electron radius (Mohr 
et al. 2012), c ^ 3.00 x 10'° cms"' is the speed of light in va- 
cuum, and Aq and fg are the transition's rest wavelength and oscil- 
lator strength, respectively; the coefficient in the equation above 
is obtained for the Nevm Aq =770.41 A transition with fo = 0.103 
(Verner et al. 1994). 

Since Tq, the equivalent width and the sensitivity in terms of 
S/N all depend on NNevm and Nncviii (see equations 1 1 and 12 below), 
our problem reduces to finding expressions for these two quantities 
in terms of and T. As a starting point, we assume that the Ne vm 
column density can be computed as 

NNcvm = I "Ncviii dl' . (2) 
Jo 

Here, / is the physical, linear extension of the gaseous, Ne viii absorb- 
ing structure (or 'cloudlet'), and wncviii is the corresponding Ne vm 
particle density, which can be expressed in terms of the Ne vm ion 
fraction, /ncviii = («Neviii/"Nc)i and the neon abundance (by number) 



' The fact that the Neviii detection limit in our synthetic spectra with 
S/N=1000 and the flattening of the Nevm CDDF both correspond to 
Nncviii ~ lO'^cm"^ is only coincidental and has therefore no deeper mean- 
ing. 
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relative to hydrogen, {nt^^/nn), as 



^Nevm 



"Nc VII 



(3) 



Let us now assume that the absorbing structures are self-gravitating 
and therefore typically have linear sizes on the order of the local 
Jeans length (Schaye 2001) 



Lj = 0.169 Mpc (;jH/10-^cm-^) 



-1/2 



X (r/10*K)'/- (/„/0.168)'''^ 



(4) 



where, fg is the gas mass fraction and fg = O^/n,,, = 0.168 its 
universal value.* 

While equation (4) is reasonable,^ we need to take into account that 
the distribution of intergalactic metals is rather inhomogeneous^ on 
scales < 10^ kpc (Schaye et al. 2003) and even on scales < 1 kpc 
(Schaye et al. 2007). We do this by introducing a factor, c/., which 
gives the average filling fraction of the metal cloudlets along the gas 
cloud, in such a way that the total intersected metal absorbing path 
length (or the integrated linear size of the metal absorber) is 



/ = ClL] . 



(5) 



Note that, by definition, < cl < 1, where cl « 1 implies a very 
inhomogeneous, and c^^ « 1 a very homogeneous metal distribution 
across the gas structure. 

If the physical conditions (i.e. temperature, density), the ion- 
isation state and the neon abundance across the metal absorber are 
uniform, then equation (2) simpUfles to 



NNevm - ikevm ( 1 Ci^Ljllji . 

Plugging equation (4) into the above equation yields 



NNevm = 5.90 X lO'^cm-^ c, iC^^™ 



(6) 



ll0-5cm-3i IsxlQSKi \ 



0.168 



1/2 



(7) 



where [Ne/H] = logjQ (iif^^/nn) - log^ (nNe/«H)o- The fiducial val- 
ues we have inserted in the above equation correspond to the peak 
ionisation fraction of Nevm in CIE, /Nevm = 0.130 (Sutherland & 
Dopita 1993), the (approximate) corresponding peak temperature, 
r = 5 X 10' K, and the solar neon abundance log,Q (nNe/nn)© = 
-3.91 given by Anders & Grevesse (1989, their Table 2). Equa- 
tion (7) readily shows that the Ne vm column density of gas with 
r * 5 X 10' K and nn ~ 10"', and a neon abundance [Ne/H] = -0.5 



is Nn 



10 cm , assuming fg - 0.168 and that the metals are 



homogeneously distributed across the structure, i.e. = 1. 

In general, both thermal and non-thermal (e.g. turbulence, 
Hubble expansion) mechanisms will contribute to the total fine 
width. If both types of mechanisms are stochastic and if they both 



* The total mass density parameter is Cl„ = Qi, + Clou- The most recent 
measurements of the baryonic and dark matter density parameters yield, 
respectively, Hj, = 0.0449 ± 0.0028 and IIdm = 0.222 ± 0.026 (Jarosik et al. 
2011). 

' We have shown (Tepper-Garcia et al. 2012; see also Altay et al. 2011; 
Rahmati et al. 2012) that our simulation (albeit at z = 0.25) is consistent 
with the assumption that typical H i absorbers are self-gravitating clouds in 
hydrostatic equilibrium with sizes ~ Lj. 

* The results by Schaye et al. (2007) correspond to observed metal-line 
absorbers atz~ 2.3. 



obey a Gaussian distribution, then the total line width may be ex- 
pressed in terms of the thermal (br) and non-thermal (bnt) broaden- 
ing as 



2kT 



+bt 



(8) 



where Boltzmann's constant k = 1.38 x 10 '^ergK and m^e = 
33.5 X 10"^ g is the mass of the neon atom (Mohr et al. 2012). The 
line broadening due to non-thermal processes, in particular turbulent 
mechanisms, might be difficult to model. However, as we show 
in Appendix B (see Figure Al), the intrinsic width of the Nevm 
absorption lines in our simulation at z = 0.5 is dominated by thermal 
broadening, i.e. br » such that the total line width may be well 
approximated by 



\l/2 



i>Nevin-Z'r = 20.3kms-'(^^^^) 
Using equations (1), (7), and (9) we get 



To(Nevin) = 3.46 cz, 10 



(9) 



■ 1 10-5cm-3 / 



0.130 
0.168 



(10) 



The above equation shows that, even for a density as low as wh ~ 
10"' cm"' and a metallicity [Ne/H] = -0.5, the opacity at the line 
centre of Ne vm absorbing gas with a temperature close to the peak 
temperature of the Ne vm ionisation fraction in CIE (r » 5 X lO' K) 
is high. To ~ 1, provided that the metals are well mixed within the 
absorbing structure (c^ ~ 1). 

Given a line strength tq, we compute the corresponding rest- 
frame equivalent width using 



%J—oo 



e-^W)rfjc, 
Toexp(-x^), X 



(11) 



(A - Af))/AA.D, and AAd 



where t(x) 
Ao{bvievm/c). 

Finally, the sensitivity in terms of the signal-to-noise ratio (S/N) 
required to detect gas at significance Ncr and the corresponding line 

strength are related by 



No- ~ A'(S/N)"' ~ (1 - e""") « To , 



To « 1 . 



(12) 



Equations (7), (9), (10), (11), and (12) thus constitute the solution 
to our problem. Of course, these equations are not restricted to 
Ne vm; they can easily be modified to describe the absorption by 
any ion of any other metal. 

In addition to the observables considered above, we provide an 
estimate of the total hydrogen column density (A'h) and the neutral 
hydrogen column density (A^hi) of the Nevm cloudlets using 



NH = nnl 
and 

A^Hi = /hi A'h 



(13) 



(14) 



where fm = (nni/nn) is the Hi ionisation fraction; this is obtained 
from the same ionisation model we use to compute the Ne vm ion- 
isation fraction (/Nevm; see Section 2). 

The actual observables, however, are the column densities 
of the gas harbouring the Nevm cloudlets, which we denote by 

= «H Lj and A'hi' = /hi N^'. Note that each of these quantities 
differs from the corresponding quantity (equations 13 and 14) by a 
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factor Cj}, i.e. A^h' = A'h c^' and Nm = Nm c^'. 

We close this section with an important remark: Strictly 
speaking, the particle density Wh and the temperature T explicitly 
entering equation (6) and equation (9), respectively, correspond to 
the particle density and to the temperature of the metal cloudlets, 
whereas the particle density and temperature implicit in equation 
(6; and all derived equations) through Lj = Ljin-a, T) refer to the 
particle density and to the temperature of the surrounding gas. 
These quantities need not be equal, in general. However, in deriving 
equations (7) and (10) we have implicitly assumed they are (at least) 
of the same order of magnitude. As we have shown previously in 
Section 2.1, this assumption is well justified. 



2.2.1 Setting the model parameters 

Our analytic model depends basically on two parameters: 1) cl, 
which measures the homogeneity of the metal distribution within 
the gas structure; and 2) the (average) neon abundance (nNe/nn) in 
the absorbing gas. The first parameter, cl, is a free parameter that 
needs to be calibrated e.g. using simulations. The second parameter, 
(rifie/nu), should be constrained by observations. Ideally, the value 
of Cl obtained from our simulation should be consistent with obser- 
vations as well; and likewise, the value of (nNe/wn) constrained by 
observations should be consistent with the typical neon abundance 
of Nevm absorbers in our simulation. 

Let us start by considering the second parameter. Figure 
7 shows the distribution of the neon abundance, [Ne/H] = 
logio (nNe/nn) - logio (mnc/wh)©. in gas in our simulation at z = 0.5, 
adopting logjg (wNe/wn)© = -3.91. The left panel shows the mass- 
weighted average of the neon abundance in gas as a function of the 
gas density and temperature. The contours in this panel indicate the 
distribution of Nevm mass in gas, and correspond to the distribution 
shown in colour in the top-left panel of Figure 2. 

The right panel displays the distribution of [Ne/H] values in 
the Nevm absorbing gas in our simulation at z = 0.5 (blue curve). 
The full range of neon abundances measured in some of the Ne vm 
absorbers observed at low redshift is indicated by the shaded area. 
The data points in both panels represent a subset of the observed 
Ne vm absorbers listed in Table 1 . This figure shows that the neon 
abundance of the Ne vm absorbers in our simulation is distributed 
over the range -1.5 < [Ne/H] < 0, with median [Ne/H] = -0.7. 
Clearly, the average neon abundance of the gas in our simulation is 
broadly consistent with the neon abundance of the observed Ne vm 
absorbers (data points) at similar densities and temperatures, al- 
though our simulation favours slightly lower values. In principle we 
could tabulate the mass-weighted mean of (nNe/nn) obtained from 
our simulation as a function of Mh and T and use it as input in our 
analytic model. However, for simplicity and as a first approach, we 
have chosen to fix the neon abundance in our analytic model to the 
median value [Ne/H] = -0.7 found in Ne vm absorbing gas in our 
simulation. This choice is supported by the fact that the metallicity 
of the warm-hot gas in our simulations is robust to changes in the 
sub-grid physics and evolves only slowly (Wiersma et al. 20 11). 

Next we need to set the value of the parameter cl. As mentioned 
before, this is a free parameter whose exact value carmot be set 
a priori. However, it is possible to constrain the range of values 
beforehand. For instance, by definition a value c « 1 implies that 
the metals are homogeneously distributed across the structure; also, 
by virtue of equations (4), (5), and (7), such a value implies that the 
sizes of absorbers with NNevm ~ 10''* are on the order of 1 Mpc, 



neither of which is consistent with observations (see Section 4). 
Thus, values in the range Cl « 1 seem more plausible. 

To determine the value of cl more precisely, we choose the 
following approach: First, we measure the Ne vm column density 
(Nnuviii), the rest-frame equivalent width (Wr) and the central optical 
depth (to) for each Ne viii absorber identified in the synthetic spec- 
tra obtained from our simulation, as explained in Appendix B. We 
use for this purpose the sample of absorbers identified in simulated 
spectra that have not been convolved with a COS-like line-spread 
function to avoid a potential bias between intrinsic/physical- and 
spectral properties (e.g. between the gas temperature and the corres- 
ponding line width; see Figure Al). Next, we extract subsets from 
our sample of Nevm absorbers chosen in terms of various NNevm, 
To, and Wr cuts, and compute for each of these Ne vm absorbers a 
characteristic (i.e. Nevm optical depth weighted) density and a 
characteristic temperature T, as described in Section 2.1. We finally 
adjust Cl in our analytic model (equations 7, 9, 10, and 1 1) to simul- 
taneously reproduce the distribution of each of these subsets on the 
na- T plane. 

We find that logiQ Cl = -0.5 (i.e. cl ~ 0.3) yields very sat- 
isfactory results, adopting our fiducial values [Ne/H] = -0.7 and 
fg = 0.168, as shown in Figure 2. The panel on the top-right of 
this figure displays a scatter plot of the Ne vm absorbers identified 
in our simulated spectra colour-coded for the various Nnc VIII cuts 
indicated in the legend. The solid contours indicate the predicted 
Ne vm column density as a function of mh and T as given by our 
analytic model (equation 7). The bottom-left and bottom-right pan- 
els display, respectively, a scatter plot of the Ne vm absorbers colour 
coded by different Ne vm central optical depth (tq) and rest-frame 
equivalent width (Wr) cuts. The contours in the corresponding panel 
show Tq as a function of nn and T as given by equation (10) and 
Wr as a function of hh and T computed using equation (11). The 
top-left panel displays the distribution of Ne vm mass in gas as a 
function of its density and temperature in our simulation and has 
already been discussed in Section 2.1. The blue curve indicates the 
minimum column density NNevm = 10"cm"^ allowed in our fitting 
procedure (see Appendix B) and has been computed using equation 
(7). Note that virtually all the Ne viii mass in gas is contained within 
the region on the n^-T plane defined by this limit. This is important 
as a consistency check because it implies that we can detect the bulk 
of the Nevm mass in absorption using our synthetic spectra. 

The contours in panels 2 - 4 of Figure 2 agree well with the 
values of the coloured points in the same panels, indicating that 
our analytic model reproduces the spectral properties of a statistical 
significant sample of simulated Ne vm absorbers as a function of 
the density and temperature of the absorbing gas. While it is true 
that our fiducial parameter values (fg - 0.168, logjQCY = -0.5, 
[Ne/H] = -0.7) may not apply in general to all absorbers, they 
should still be representative. First, there is no particular reason why 
the ratio of baryonic to total mass in such systems should be very 
different from its universal value. Moreover, the value logjQ Cl = 
-0.5 yields sizes for the Ne vm absorbers (see equation 5) on the 
order of 10 - 10^ kpc, which are consistent with the typical size 
of Ne vm absorbers around haloes both in observations (Meiring 
et al. 2012) and in other simulations (Ford et al. 2012). Values Cl ^ 
1 are also consistent with observations of C iv absorbers at high 
redshift (Schaye et al. 2007). And finally, the fiducial value of the 
neon abundance is broadly consistent with the currendy measured 
abundances of Ne vm absorbers at low redshift. 

In principle, higher values of [Ne/H] and lower values of Cl 
are both allowed by observations. Within the framework of our 
model, however, these parameters are not completely independent 
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Figure 7. Neon abundance, [Ne/H] = logm (hnc/ih) - logio ("Ne/^n)©. in g^** in our simulation at z = 0.5. Left: Mass-weighted mean of [Ne/H] in gas as a 
function of tlie gas density and gas temperature. The black contours indicate the distiibution of Ne viii mass in gas, and they correspond to the distribution 
shown in colour in Figure 8. Right: Distribution of [Ne/H] values in Ne viii absorbing gas in our simulation (blue curve). The vertical dashed line indicates the 
corresponding median value [Ne/H] = -0.7, which is also the fiducial value adopted in our analytic model. The shaded area indicates the full range of neon 
abundances measured in Ne viii absorbers at low redshift (see Table 1). The data points in both panels represent a subset of the observed Ne viii absorbers listed 
in Table 1 . Note that their measured neon abundances (right panel) are broadly consistent with the mean neon abundance of the gas in our simulation at the 
corresponding densities and temperatures (left panel). All abundances are given relative to logjo («Ne/nH)o = -3.91. 



of each other. For instance, if we demand that our analytic model 
yields quantitatively the same results in order to remain consistent 
with our simulation, then a higher neon abundance requires a lower 
value of cl for a given ;?h and T. In other words, our analytic model 
predicts that at a given Ne vin column density, high (low) metallicity 
absorbers are more compact (extended), on average. 

Needless to say, it is not the exact value of Cl (nor of [Ne/H]), 
but its order of magnitude that is important as a parameter that 
quantifies the degree of metal mixing in gas. The power of our 
analytic model ultimately lies in providing a connection between the 
observable parameters (Nncviii- ^r^Nm) and the physical parameters 
(wh, T, Nh, [M/H], /) of the absorbing gas, while accounting for the 
observed inhomogeneous distribution of metals. 



3 SPECTRAL SIGNATURES OF Ne viii ABSORBING GAS 

Now that we have found sensible values for the parameters entering 
our analytic model, we use it, in combination with our simulation, 
to analyse the connection between the spectral characteristics and 
the physical state of Ne vm absorbing gas at z = 0.5. The results are 
shown in Figure 8 and discussed below. 

Each of the panels in this figure displays a series of contours 
that show the behaviour predicted by our analytic model for Ne vin 
absorbers as a function of the density and temperature of the Ne vin 
absorbing gas. The physical parameter we consider are: the Ne vin 
column density (top-left), the Ne viii rest-frame equivalent width 
(top-right), the Ne vm central optical depth (middle-left), the min- 
imum S/N (at 5o" significance) required to detect Ne vm absorbing 
gas assuming the line to be fully resolved (middle-right), the total 
hydrogen column density Nh (bottom-left), and the neutral hydro- 
gen column density A'hi (bottom-right) of the Nevm bearing gas. 
The alternative y-axis in the top-right panel indicates the thermal 
broadening of a Nevm absorption line (equation 9). The dashed 
lines in the middle-right panel indicate the approximate size (in kpc) 
across the Ne vm absorber, given by equation (5). The column dens- 
ity, the central optical depth, and the rest-frame equivalent width are 
computed using equations (7), (10), and (11), respectively, adopting 



logio ct = -0.5 (i.e. Cl ~ 0.3) and the fiducial values [Ne/H] = -0.7 
and fg = 0.168. The coloured areas (which are identical in all pan- 
els) show the distribution of Ne vm mass in gas in our simulation at 
z = 0.5 and have already been discussed in Section 2.1. Note that 
they are identical to the distribution shown in colour in the top-left 
panel of Figure 2. 

We first give a few comments on the qualitative trends seen 
in each of the panels in Figure 8. Quite apparent is the similarity 
between the behaviour of the Ne vm ion fraction (/ncviii; Figure 1) 
and the behaviour of each of the line observables with iih and T, 
indicating that /Nevm is the dominant factor in setting the absorption 
signatures of the absorbing gas. Yet, the correspondence is not 
perfect. In particular, we see that despite the relatively high ion 
fraction /ncviii ~ 0.1 of low-density, photo-ionised gas with «h ~ 
IO"*cm"^ and T ~ W* K, the corresponding Ne vm column density 
is low, Nncviii ~ lO'^cm"^, compared to the Nevm column density 
NNeviii ~ 10'''cm-2 of gas with iia ~ lO^^cm-^ and T > 10^ K. This 
cab be understood in terms of equation (7), which predicts that, for a 
fixed /ncviii and a fixed neon abundance, the Ne vm column density 
scales as Nncviii «^ riyi I oc (im T)''^. As a result, the Nevm column 
density of photo-ionised gas at «h < 10"^ cm"' and T < 10^ K can 
be orders of magnitude lower than the Ne vm column density of 
collisionally ionised gas with > lO^^cm"' and T x 5 x 10^ K, 
even though /ncviii ~ 0. 1 in both gas phases and even if we assume 
identical metallicities. The equivalent width (IV, ; top-right panel) of 
the Ne vm absorbing features produced by these gas phases behaves 
in a similar way, since for NNevm S 10'*cm"^, which is the case of 
much of the gas in the range of densities and temperatures shown, 
Wr is directly proportional to NNevm (see right panel of Figure Al). 

The optical depth (middle-left panel) of the gas shows a slightly 
different behaviour. For instance, the optical depth at the line centre 
of Ne vm absorbers produced in photo-ionised and collisionally ion- 
ised gas is comparable, Tq > 0.1, even though the corresponding 
column densities are an order of magnitude apart. This can be ex- 
plained as follows: At temperatures T < 10' K the Ne vm absorption 
lines are expected to be narrow (iiNcvm ^10 km s"' ; see right y-Rxh 
in top-right panel) if their width is dominated by thermal broadening. 
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Figure 8. Ne viii line observables and mass distribution of Ne viii in gas at z = 0.5 as a function of hydrogen particle density (jih) and temperature (T). The 
coloured regions show the fraction of Ne vra mass in the simulation indicated by the colour bar on the bottom-right margin, and they are identical in all panels 
(and correspond to the coloured regions shown in Figure 2). The contours in each panel show the predictions of our analytic model. Top-left: Ne vra column 
density (Nnc viii in cm"^) as given by equation (7). Top-right: Rest-frame equivalent width (in mA) of a Ne vm absorption line computed using equation (11). 
The right i/-axis in this panel indicates the corresponding thermal broadening (equation 9). Middle-left: Ne viii central optical depth (tq) as given by equation 
(10). Middle-right: Minimum S/N required to detect Ne vin absorbing gas calculated assuming a 5o" significance, assuming the line to be fully resolved. The 
dashed contours indicate the approximate size (in kpc) across the Ne vin absorber. Bottom-left: Total hydrogen column density (A^h in cm"^) of Ne vm bearing 
gas (equation 13). Bottom-right: Neutral hydrogen column density (Wm in cm"^) of Ne viii bearing gas (equation 14). Note that the total H and Hi column 
densities of the gas clouds that harbour the metal patches ai'e a factor c^' greater. The contours shown in all panels have been computed assuming a neon 
abundance [Ne/H] = -0.7 (relative to logjQ (nNe/'Hl)© = ~3.91), a gas fraction fg = 0.168, and a line of sight metal filling fraction logmf/. = -0.5 (i.e. 
ct « 0.3). Note that the contour values are given on a logarithmic scale in all the panels on the left and in the bottom, and on a linear scale in the others. 
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Since To oc (N^evm/^Neviii). a decrease in Nncviii is counterbalanced 
by a decrease in Z^Nevm such that tq remains approximately constant. 
Since the minimum S/N value required to detect Ne vin in gas at 
a given density and temperature scales as (tq) ' (see equation 12), 
there is a natural correspondence between low (high) optical depths 
of the Ne vm absorbing gas and the high (low) S/N values required 
to detect such gas, as can be seen by comparing the middle panels 
with each other. 

Finally, the total hydrogen column density A^h scales as 
(whT)''^ (see equation 13), which explains the linear (increasing) be- 
haviour of logio A'h with logio ita (logm T) for a fixed T (nn) as seen 
in the bottom-left panel. Likewise, the typical Unear size across the 
Ne vm absorber shown by the blue dashed lines in the middle-right 
panel obeys I oc Lj oc (T/nuY^^, thus linearly increasing (decreas- 
ing) with logiQ T (logiQ^H)- The behaviour of N^i is dictated by 
the non-trivial dependence of /hi on hh, and particularly on T (see 
equation A7 of Schaye 2001, and Figure 7 of Tepper-Garcia et al. 
2012). 

Further insight into the nature of Ne vm absorbers at low red- 
shift can be obtained by comparing the distribution of Ne vm mass in 
gas in our simulation (coloured regions) with the spectral properties 
of the gas as a fimction of its density and temperature predicted 
by our analytic model (contours) shown in each panel. We see that 
in our model universe most of the Ne vin mass is contained in gas 
with densities logiQ(«H/cm'') < -4 (corresponding to overdensit- 
ies A < 150 at z = 0.5) and temperatures T < 10* K for which 
the Ne VIII column densities are expected to be N>jeviii ^ lO'^^cm^^ 
and the corresponding rest-frame equivalent widths Wr < 30 mA. 
This result is consistent with the distribution of column densities 
and equivalent widths obtained from our synthetic spectra (Section 
2.1.2). The gas containing most of the Ne vm mass in our simulation 
thus produces weak to moderate Ne vm absorption features with 
optical depths at the centre log|Q To < -0.5; therefore, S/N values of 
30 or higher are required to detect these features (at 5cr significance). 
The Une widths implied by the corresponding gas temperatures, if 
thermal broadening dominates, are ^Nevm < 30km s~\ with most 
of the Ne vm mass in gas producing thermally broadened lines with 
lOkms"' < ^Nevm < 30km s '. The gas responsible for the Nevm 
absorption produces Hi absorption displaying a wide range of Hi 
colimin densities: 10" < (iVni/cm"^) < 10'**. The average baryon 
content of this gas in the form of hydrogen corresponds to total 
column densities log|Q(A'H/cm"^) < 19.5. 

We find that only a small fraction of the Nevm mass in 
our simulation resides in gas with densities n-a ^ 10"''cm"^ and 
temperatures around T 5 x 10' K, i.e. in gas which produces 
Nevm absorption with column densities NNevm ^ lO'^cm"^ (and 
Wr > 30 mA). These Ne vm absorbers should be detectable in ab- 
sorption in spectra with S/N > 10 (at 5cr significance), given that 
their line strength is log,,, Tq > -0.5. The gas temperatures imply 
line widths foNevm ^ 20 km s"' . The average total hydrogen column 
densities of the gas responsible for this type of Ne vm absorption 
can be as high as A^h ~ lO^^cm"^. Due to the high gas temperat- 
ures, the hydrogen neutral fraction in this gas is low, and the gas 
produces H i absorption of 'BLA' type: weak, with column densities 
10" < (WHi/cm-2) < 10■^ and broad, withi^Hi i 90 km s"' ; a high 
sensitivity is thus required to detect this gas phase using H i absorp- 
tion. The linear extent of absorbers with densities itn ~ lO'^'cm"' 
and temperatures J < 10'' K is on the order of 10^ kpc, and it de- 
creases with density. Note, however, that the size, and the H and H i 
column densities of the gas clouds that harbour the metal cloudlets 
responsible for the Ne vm absorption are a factor c^' higher. It is 
worth noting that all these properties closely resemble the properties 



of the Ne vm absorbers at low redshift reported in the literature to 
date. We will discuss this in more detail in the next section. 

Our analytic model, in combination with our simulation, there- 
fore predicts that Nevm absorbers with NNevm ^ 10''*cm"^ (and 
Wr ^ 30 mA) are not common at low redshift. However, the weakest 
members in this class of Ne vm absorbers are just strong enough to 
allow their detection in absorption spectra with S/N ~ 10. The detec- 
tion of weaker, but more common Ne vm absorbers requires a signi- 
ficantly higher sensitivity S/N ~ 100. Our model thus simultaneously 
provides a natural explanation for the low detection rate of Ne vm 
absorbers at low redshift -first noted by Narayanan et al. (2012)-, 
as well as for the similarity in their strength (Nncviii ~ 10 ''^cm"^) 
given that the typical sensitivity of common UV spectra is S/N ~ 10 
(see Section 4). 

Equally important, the result from our model that absorbers 
with NNevm ^ lO'^^cm"^ and Wr > 30 mA are only produced in gas 
with temperatures in a very narrow range around T « 5 x 10' K 
strongly supports the idea that this particular class of Ne vm ab- 
sorbers is a sensitive probe of the gas temperature (Savage et al. 
2005). In other words, the detection of Nevm absorbers with 
NNevm ^ lO'^'cm"^ is a sttong indication of the presence of gas 
with temperatures logjQ(r/ K) ^ 5.5, under ionisation equilibrium 
conditions. 

Our model further predicts that the gas that produces such 
strong Ne vm absorption can have total hydrogen column densities 
A'h » 10"cm"^. Nevertheless, the overall baryon content of low- 
redshift Ne vm absorbers may not be as high as had been expected 
(e.g. Savage et al. 2005). Indeed, the total baryon content of the 
gas producing the Ne vm absorption in our simulation amounts to 
D;XNe viii) = 1.22 x 10"^ or roughly 3 per cent of the total baryon 
budget. The gas producing absorbers with NNevm ^ 10'"* cm"^ con- 
tain less than 0.5 per cent of the total baryons. Note, however, that 
our estimate represents a strict lower limit, since we ignore absorbers 
with Nnctoi < lO'^cm"^ which contain roughly 20 per cent of the 
Nevm mass in our simulation (see Section 2.1.2), although accord- 
ing to our model predictions their average total hydrogen column 
densities, and in consequence their baryon content, is expected to 
be relatively low, A'h < 10"*cm"^ . Note also that the total baryon 
density traced by the Ne vm absorbers is a factor c^' higher, or ap- 
proximately 10 per cent of the total baryon budget in our simulation 
if Cl « 0.3. Our result is of the same order as the corresponding 
value nj(Ne vm) = 2.63 x 10"^ (or 6 per cent of the total baryons in 
the Universe if fl/, = 0.0418) estimated by Narayanan et al. (2009), 
based on a single detection, and the value ^^(Ne vm) = 2.00 x 10"' 
(roughly 4 per cent of the baryon density if f2j, = 0.0418) obtained 
by Meiring et al. (2012) based on three detections. 

On the other hand, if < 0.3, then the baryon content traced by 
Ne vm in our simulation could be significantiy higher. For instance, 
if cl ~ 0.1 as is suggested by observations in combination with our 
analytic model (see Section 4), then the fraction of the total baryon 
contained in the (warm-hot) gas harbouring the Ne vm would be 
boosted to roughly 30 per cent, which is comparable to the fraction 
of baryons in gas traced by BLAs in our simulation (25 per cent; 
Tepper-Garcia et al. 2012). In this case, however, it is unclear what 
the overlap between BLAs and Ne vm is in terms of the baryons 
traced by these ions, and in consequence the true baryon content of 
warm-hot diffuse gas at low redshift remains uncertain. 
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Table 1. Summary of Ne viii detections at z < 1. See text for details on the notes. List of references: (1) Savage et al. (2005); (2) Narayanan et al. (2009); (3) 
Narayanan et al. (2011); (4) Narayanan et al. (2012); (5) Meiring et al. (2012); (6) Savage et al. (2011); (7) Tripp et al. (2011). 
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[cm-2] 


[kms"'] 


[mA] 


[cm-3] 


[K] 




[cm-2] 


[cm-2] 


[kpc] 




HE 0226-4110 


0.207 


13.89 ±0.11 


22.6 ± 15.0 


36 ±10 


-3.20 


5.68 


-0.25 


20.06 


13.87 


60 


(1,6) 


3C 263 


0.326 


14.06 ± 0.08 


52.4 ± 8.8 


57± 10 


-2.79 


5.72 


0.04 


19.48 


13.09 


6 


(2,4) 


PKS 0405-123 


0.495 


13.96 ± 0.06* 




45 ±6* 


-3.10'' 


5.67'' 


-0.44 


19.70'' 


13.50 


20^ 


(3) 


PG 1148+549 


0.684 


13.95 ± 0.04 


32.0 ±5.0 


43 ±4 


-3.09 


5.68 


-0.36 


19.80 


13.65 


10 


(5) 


PG 1148+549 


0.705 


13.86 ± 0.05 


27.8 ±3.5 


35 ±4 


-2.34 


5.69 


0.14 


19.15 


12.90 


1 


(5) 


PG 1148+549 


0.725 


13.81 ±0.07 


41.4 ±7.5 


33 ±5 


-3.96 


5.72 


0.14 


18.83 


13.08 


20 


(5) 


PG 1206+459 


0.927 


13.71 ±0.29 








5.56 




19.80 






(7) 


PG 1206+459 


0.927 


14.04 ± 0.08 








5.60 




20.30 






(7) 


PG 1206+459 


0.927 


14.53 ± 0.04'' 








5.65 




20.50 






(7) 



" Measured using Gaussian profile fitting, unless stated otherwise. 

* Measured using the apparent optical depth method (AOD; Savage & Sembach 1991). 

Computed here using the measured column density and Doppler parameter assuming a Voigt-profile, unless stated otherwise. 
'' Estimated assuming photo- and collisional ionisation equilibrium. 
^ Inferred from the total hydrogen column density and nn using I = Nu/ir. 
f Estimated here from the inferred size / of the photo-ionised phase and Nyi using tin = Nji/l- 
^ Relative to logio ("Ne/in)© = -3.91. 

* Total column density of a two-component feature. 



4 COMPARISON WITH OBSERVATIONS 

Given the relevance of the results presented in the previous section, it 
is important to test the predictions from our analytic model using the 
measured (and inferred) properties of Ne viii absorbers observed at 
low redshift. In Table 1 we present all the Ne vm detections available 
in the Uterature to date, listing only those physical parameters that 
are relevant for our study (if available). Unless stated otherwise, 
all the values listed in Table 1 are taken from the corresponding 
reference given in the last column. 

In some cases (e.g. Savage et al. 2005; Meiring et al. 2012), the 
gas temperature of the coUisionally ionised gas phase responsible 
for the Ne vm and O vi absorption has been obtained by scaling A'h 
in order for the model to reproduce the measured Ne vm and O vi 
colunm densities, assuming CIE and some value for the metallicity 
(usually the same value of the metallicity inferred for the photo- 
ionised gas phase). In these cases, A'hi has been calculated from 
the inferred temperature which sets the Hi ionisation fraction /hi 
(assuming CIE) and the inferred A'h using A^hi = /hi A^h- 

In others cases (e.g. Narayanan et al. 2012), the gas temperat- 
ure of the observed Ne vm absorbing gas has been inferred from the 
ratio of the observed Ne vm to O vi column densities (NNevm/Novi) 
assuming CIE and fixed (e.g. solar) abundances. The inferred tem- 
perature (together with the assumption of CIE) has then been used 
to convert the measured neutral hydrogen colunm density (A^hi) into 
a total hydrogen column density (A'h). In these cases, the H i column 
density of the gas phase responsible for the Ne vm and O vi absorp- 
tion has been either directly measured in the spectrum (e.g. Savage 
et al. 2011) or estimated by overlaying onto the Hi absorption fea- 
ture detected in the observed spectrum an additional H i component 
whose width is set by the inferred gas temperature and whose H i 
colunm density is the largest possible such that the overall fit is still 
consistent with the data (e.g. Narayanan et al. 2011, 2012). 

With one notable exception, there are no estimates (apart 
from upper or lower limits) of the density of the Ne vm absorbers 
available in the Uterature. The only estimate is given by Naray- 



anan et al. (2011) for a Nevm absorber at z ~ 0.5. They obtain 
logio(nH/cm""') = -3.10 using a model that combines the effects of 
photo- and collisional ionisation (as opposed to assuming CIE only) 
to compute the ionisation balance of the absorbing gas, which is the 
same approach we use here (see Appendix A for more details on the 
analysis of this absorbing system). We therefore consider their meas- 
urements to be the most suitable for comparison with the results 
from our analytic model. As for the other observed Ne viii absorbers, 
and with no other choice at hand, we compute the gas density from 
the inferred size (/) and the total hydrogen density (A^h) using the 
simple relation mh = Nn/l, whenever possible. Note that I refers 
to the linear extension of the photo-ionised component inferred by 
assuming photo-ionisation equilibrium (PIE) while A'h refers to the 
total hydrogen column density of the associated coUisionally ionised 
gas phase. 

Finally, in order to be self-consistent, we compute the rest- 
frame equivalent width (rather than using its measured value) with 
the help of equations (1) and (11) using the values for the Nevm 
column density (Nncviii) and the Doppler parameter (Z^Nevm) ob- 
tained by each group from line-profile fitting. 

Figure 9 shows the comparison of the physical parameters 
predicted by our analytic model to the corresponding measured and 
inferred parameters of the Nevm absorbers listed in Table 1. Note 
that we use only the subset of the observations for which we have 
an estimate of the gas density (nn)- The panels in this figure display 
a series of contours which show, from the top-left and moving 
clockwise, the Ne vm column density (NNevm), the rest-frame Ne vm 
equivalent width (Wr), the total hydrogen column density (A'h; black 
lines) and the neutral hydrogen column density (A'hi; blue curves) 
of the Ne vm absorbing gas, and the minimum required sensitivity 
(in terms of S/N) to detect gas using Ne vm absorption as a function 
of the gas density and temperature, assuming the line to be fully 
resolved. The top-right panel shows additionally the approximate 
size (in kpc) across the Nevm absorber (blue dashed lines). Note 
that the contours are identical to the corresponding contours shown 
in Figure 8. The data points represent the observed Ne vin absorbers. 
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Figure 9. Comparison of the spectral and physical properties predicted by our analytic model to the corresponding measured and inferred properties of a sample 
of observed Ne vm absorbers at low redshift as reported in the literature: Nevm column density (NNeviii; top-left), Nevni rest-frame equivalent width (W,-; 
top-right; black contours) and linear size in kpc (/; top-right; blue dashed lines), total hydrogen column density (A'h; bottom-right; black lines) and neutral 
hydrogen column density (Wm; bottom-right; blue curves). The contours in the bottom-left panel show the minimum sensitivity in terms of S/N required to 
detect gas at a given density and temperature using resolved Nevm absorption at Scr significance. The values for each quantity measured (or inferred) from 
observations are indicated in the legend along with their coiTesponding reference. A complete list of the quantities measured and infen'ed from observations 
is given in Table 1 . The S/N values quoted in the bottom-left panel correspond to the average S/N of the data in the spectral range of interest. Note that the 
contours in each panel are identical to the corresponding contours shown in Figure 8. List of references; (1) Savage et al. (2005); (2) Narayanan et al. (2009); (3) 
Narayanan et al. (201 1); (4) Narayanan et al. (2012); (5) Meiring et al. (2012); (6) Savage et al. (201 1). 



In each panel, the measured (or inferred) value of a given physical 
parameter (e.g. NNevm) is indicated in the legend, along with the 
reference it has been taken from (see the figure caption and Table 
1). 

We find broad agreement between the results inferred from ob- 
servations and the predictions from our analytic model. The panels 
in the top row, for instance, demonstrate that the measured Ne vin 
column densities, as well as the measured rest-frame Ne vin equival- 
ent widths, of the observed Ne viii absorbers with a given inferred 
density and inferred temperature fall within the range of Ne vin 
column densities and equivalent widths predicted by our model 
for gas in the corresponding density and temperature range. It is 
remarkable that all the observed Ne vin absorbers have very similar 
strengths, NNeviii ~ 10'''cm"^, within the given uncertainties, and 
that they have all been detected in spectra with (average) S/N ~ 10, 



which is perfectly consistent with the result from our model (com- 
pare top-left and bottom-left panels). The bottom-right panel shows 
that both the ranges of H and H i column densities of the gas respons- 
ible for the Ne vm absorption obtained from observations are well 
predicted by our analytic model. Also, the typical linear extension 
of the absorbers resulting from our model is broadly consistent with 
the sizes inferred from observations (see top-right panel). The worst 
agreement between our model predictions and the corresponding 
results obtained from observations, in particular for Z, is shown by 
the Ne vm absorbers at z = 0.705 and z = 0.725 reported by Meir- 
ing et al. (2012, left and right red dots). Some of the discrepancy 
may come from the fact that our model curves are computed assum- 
ing [Ne/H] = -0.7, whereas these two observed Nevm absorbers 
have higher (super-solar) abundances [Ne/H] = 0.14 (relative to 
logiQ (nNe/^n)© = -3.91). In contrast, all the measured properties 
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of the z = 0.495 Nc viii absorber detected by Narayanan ct al. (201 1, 
black square) are almost perfectly reproduced by our analytic model, 
although the neon abundance in this absorbers is also inferred to be 
higher, [Ne/H] = -0.44, than our adopted value, [Ne/H] = -0.7. 

The fact that the average neon abundance of the observed 
Neviii absorbers ([Ne/H] ~ -0.2; see Figure 7) is higher by 0.5 
dex than the fiducial value of our analytic model has an interesting 
imphcation. As discussed in Section 2.2.1, consistency between our 
analytic model and our simulation demands that cl decreases (in- 
creases) if [Ne/H] increases (decreases). More precisely, increasing 
the value of [Ne/H] by 0.5 dex higher would require us to decrease 
logiQ Cl by the same amount. Therefore, adopting [Ne/H] = -0.2 
would require log^ Ci = -1. In other words, the typical neon abund- 
ance observed in low redshift Ne viii absorbers, in combination with 
our analytic model, favours values ci ~ 10"', thus supporting the 
idea that metals are in-homogeneously distributed at low redshift. 
Such a change in the values of [Ne/H] and Cl would leave the 
predicted values for NNevm and (top-row panels of Figure 9) 
unaltered, but it would reduce the predicted sizes /, and the gas 
column densities A'h and A'hi by a factor 10"^ » 3, all of which 
would still be consistent with the results inferred from observa- 
tions. Note that, in particular, the predicted sizes of absorbers with 
NNevm ~ 10'*cm"^ would be on the of order 10 kpc. 

It is noteworthy that the inferred gas temperature of all observed 
Ne vm absorbers is very similar and close to the temperature T « 
5 X 10' K at which the Nevm fraction peaks in CIE. In contrast, 
their densities span almost two orders of magnitude. However, if 
we ignore the metal-rich Ne vin absorbers detected by Meiring et al. 
(2012), then all Nevm absorbers observed at low redshift appear 
to be produced under rather uniform physical conditions, that is in 
gas with rin ~ 10"^ cm"^ and T « 5 x 10' K, and total hydrogen 
column densities A^h ~ lO^^cm"^. The relatively high densities of 
the observed Ne vm absorbers, which correspond to overdensities 
of A K 1.5 X 10^ at z = 0.5, indicate that they likely originate in the 
immediate vicinity of galaxies. 

It is important to stress that we do not expect a perfect agree- 
ment between our analytic model and the observations for several 
reasons. First, an accurate comparison of the predicted and meas- 
ured column densities and equivalent widths is difHcult, given that 
our model assumes the lines to be fully resolved while the intrinsic 
velocity structure of the observed lines may be blurred by the lim- 
ited resolution and sensitivity (S/N) of the data, as demonstrated 
for instance in Figure 4. Second, the physical parameters of the 
Ne vm absorbers obtained from observations are based on models 
that assume collisional ionisation equilibrium (but see Narayanan 
et al. 201 1), while we take both the effect both of photo- and colli- 
sional ionisation into account to compute the ionisation balance of 
the absorbing gas (although we also assume ionisation equilibrivun). 
Furthermore, the gas densities (ne) of the observed Ne vm absorbers, 
which we have computed from the inferred total hydrogen column 
densities (A'h) and the size (/) of the photo-ionised gas phase, are 
uncertain due to the uncertainties in both A'h and /. Also, the sizes 
predicted by otu" model correspond to those of the metal cloudlets, 
and these may be smaller than those of the gas clouds in which the 
metal cloudlets are embedded if the metals are poorly mixed on 
small scales, as we have assumed here. Finally, our analytic model 
assumes local hydrostatic equilibrium (Schaye 2001) which may 
hold only approximately and not necessarily for all absorbers. 

In spite of these uncertainties, we find that our analytic model 
broadly reproduces the properties of observed Ne vm absorbers at 
low redshift, although the number of current detections is still small. 
Thus, more data is required before we can test the predictions of 



our model in a more robust way. Such data should soon become 
available from the analysis of a large sample of high-quality spectra 
obtained with COS. 

The results presented above lead us to conclude that the Ne vm 
detections reported to date in the literature likely correspond to colli- 
sionally ionised, baryon-rich (in terms of A^h with A'h ~ 10^**cm"^), 
compact (~ 10 kpc) metal absorbers produced in gas with temper- 
atures well above 10' K and with densities mh s> 10^'cm"^. These 
strong Nevm absorbers with NNevm ~ lO'^'cm"^ most probably ori- 
ginate close to galaxies, and they make up only a small fraction of 
the general Ne vm absorber population at low redshift. The latter 
consistsmainly of weaker Nevm absorbers with NNevm < lO'^^cm"^, 
which are produced in more diffuse gas with hydrogen column dens- 
ities A'h ^ 10"cm"^, and whose detection requires a sensitivity in 
terms of S/N on the order of 100. 



5 SUMMARY AND CONCLUSIONS 

In recent years a substantial fraction of the ultra-violet (UV) spec- 
troscopy of quasars (QSO) has been dedicated to tracking down 
a peculiar gas phase characterised by its low densities (nn ~ 
10"* - lO^'cm"'') and relatively high temperatures (T > 10' K) 
which is expected to harboiu' a substantial fraction of the baryons 
in the Universe at low redshift. Next to other ions such as five- 
times ionised oxygen (O vi) and broad H i Lyo- absorbers (BLAs), 
seven-times ionised neon (Ne vm) has been considered to be a po- 
tential tracer of this gas phase, especially because of the temperature 
(r « 5 X 10' K) at which the Nevm ionisation fraction is highest 
in collisional ionisation equilibrium. But despite the relatively high 
abundance of neon in the Universe, there have so far only been a 
small number of Ne vm detections at low redshift, and consequently 
the physical conditions and the fraction of the baryons traced by 
the gas phase producing the observed Ne vm absorption are still 
uncertain. 

In this paper we have investigated the physical state of the 
Nevm absorbing gas at low redshift following two independent 
but complementary approaches. First, we have used a large cosmo- 
logical simulation of structure formation which includes many of 
the physical processes relevant to the production and distribution 
of metals in the Universe. Second, we have developed an analytic 
model to predict the signatures of Ne vm absorbing gas as a function 
of the gas density, temperature, and the abundance of metals in 
the gas. The assumptions underlying our model are that the metals 
are in-homogeneously distributed within the gaseous structures re- 
sponsible for the absorption, that these structures are self-gravitating 
with sizes of the order of the local Jeans length, that the gas is only 
exposed to the UV/X-ray background, and that it is in ionisation 
equilibrium. The last two assumptions were also used in the simula- 
tion. We have shown that the simulation and the analytic model are 
consistent with each other, and that the predictions from our analytic 
model are broadly consistent with the results from observations, 
provided that the parameter quantifying the degree of metal mixing 
along the sightline across the structure, ci^ ~ 10"'. The parameter 
Cl corresponds to the fraction of the length of the gas cloud along 
the line of sight that intersects cloudlets enriched with neon, where 
the enriched cloudlets have a similar density and temperature as the 
surrounding gas. 

The most important results we obtain from the combination of 
our simulation and our analytic model regarding the physical state 
of the Ne vm absorbing gas at redshift z = 0.5 are the following: 
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(i) High Ne vm ion fractions /n^ vm -0.1 can occur both in photo- 
ionised and collisionally ionised gas, but most of the Ne vm in our 
simulation is produced by coUisional ionisation in gas with densities 
nn ~ lO-'cm-' and T « 5 X lO^ K (Section 2.1). 

(ii) Nevin absorption with Nevm column densities NNevm ^ 
lO'^'cm"- (or Wr > 30 mA) is only produced in collisionally ionised 
gas with temperatures around T = 5 x 10^ K and Wh ^ 10^''cm~^, 
total hydrogen column densities A?h ~ lO^^cm"^, and neutral hydro- 
gen column densities A^hi ~ 10" - lO'^'cm"^. The detection of such 
absorbers is therefore a strong indication of the presence of gas with 
T > 10' K, provided the gas is in ionisation equilibrium (Section 3). 

(iii) At z = 0.5 the number of strong Ne vm (WV > 30 mA which 
corresponds to NNevm ^ lO'^'cm"^) absorbers per unit redshift is 
predicted to be dN/dz ~ 10"', which is low compared to the total 
number density of Ne vm absorbers in our simulation, dN/dz(Wr > 
0.3 mA) ~ 4 1 . The more general population of Ne vm absorbers at 
low redshift consists of weaker absorbers with NNevm 10'''cm"^ 
whose detection requires S/N > ICQ (Sections 2.1.2, 3). 

(iv) At a given Ne vm column density, the linear extent of the ab- 
sorbers scales inversely proportional to the metal abundance. Metal- 
rich ([Ne/H] « 0), strong (NKevm ^ lO'^'cm"-) absorbers have 
typical sizes ~ 10 kpc (Sections 2.1.1, 2.2.1, 4). The gas clouds 
harbouring these Nevm cloudlets are, however, a factor c^' larger. 

(v) The gas producing the Ne vm absorption in our simulations 
contains roughly 3 per cent of the total amount of baryons; absorbers 
with NNevm ^ 10"cm"^ contain less than 0.5 per cent of the total 
baryons in our simulation. The baryon content of the gas harbouring 
the Ne vm absorbers (and which has a similar density and temper- 
ature as the Nevm cloudlets) is in each case a factor c£' higher 
(Section 3). 

(vi) The neon abundances of the Nevm absorbing gas in our 

simulation are distributed in the range -1.5 ^ [Ne/H] < 0, with 
median [Ne/H] = -0.7 (Section 2.2.1). 

(vii) The distribution of Ne vm column densities obtEiined from 
synthetic spectra at z = 0.5 is well described by a single power 
law with index j8 = 2.9 ± 0.1, in the range logiQ(NNevm/cm"^) £ 
[12.5,14.5] (Section 2.1.2). 

(viii) The small number of Ne vm absorbers observed at low 
redshift to date appear to be part of a particular class of scarce, 
strong (NNevm ~ lO^^'cm"^), baryon-rich (A'h ~ lO^^cm"^), compact 
(~ 10 kpc) metal absorbers produced in gas with densities wh ~ 
10~^cm"^ and temperatures T « 5 x 10' K. The relatively high 
densities of these Ne vm absorbers indicate that they likely originate 
in the immediate vicinity of galaxies (Section 4). 

We conclude that strong Ne vm absorbers are robust probes of 
shock-heated diffuse gas. Spectra with signal-to-noise ratio ~ 10^ 
would allow the detection of the weaker systems that trace a sub- 
stantial fraction of the baryons at low redshift. 
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APPENDIX A: CURRENT STATUS OF Ne vin 
OBSERVATIONS AT LOW REDSHIFT 

In this section we give a summary of Xhe. intervening Neviii 
absorbers detected at low redsliift and reported in the literature to 
date. So-called 'proximate' Ne vin absorbers, which are believed to 
originate in the outflows of QSOs, have been also identified at low 
redshift (e.g. Muzahid et al. 2012), but they are out of the scope of 
this study. 

The first Nevm absorber (with a total column density 
log]Q(NNeviii/cm"^) = 13.85 and corresponding rest-frame equival- 
ent width Wr ~ 33 mA) was discovered by Savage et al. (2005) 
at z = 0.207 along the sightline to the QSO HE 0226-4110 

= 0.495) in a combined FUSE / STIS spectrum with S/N 
= 10-20 per resolution element. Absorption by H i and low ionisa- 
tion states of other elements (e.g. C m, N m. Si iii) and highly ionised 
oxygen (O vi) where also detected at the same redshift. While the 
low ionisation states could be explained by photo-ionisation, the 
same model failed to explain the observed Ovi - Nevm absorp- 
tion. The implied path length (~ 10 Mpc), in particular, would 
result in a Hubble broadening which is 10 times larger than the 
observed line width. Savage et al. (2005) thus concluded that the 
O VI - Ne vm absorption is produced by coUisionally ionised gas at 
a temperature logjo(r/ K) « 5.7. They estimated a total hydrogen 
column density log|Q(7VH/cm"^) ~ 20 for this gas phase, suggest- 
ing that Nevm absorbers may harbour a significant fraction of the 
baryons in the Universe at low redshift. This absorption system 
was re-observed by Savage et al. (201 1) in a COS spectrum of 
the source HE 0226-4110 with S/N=20 - 40. With the Nevm ab- 
sorption out of the COS spectral range at this redshift, this group 
focused on the detection of a feature identified as a broad Lya 
absorber (BLA) with logio(A^Hi/cm~^) « 13.9 and a Doppler para- 
meter i>Neviii = 72^" kms"', consistent with gas at a temperature 
in the range logio(r/ K) = 5.4 - 5.7. A joint model of the BLA - 
O VI - Ne VIII system allowed them to constrain the temperature to 
logiQ(r/ K) = 5.68 ± 0.02 and the total hydrogen column density to 
logiQ(A'H/cm"^) = 20.06 ± 0.09, thus reinforcing the conclusion by 
Savage et al. (2005) that Ne viii absorbers at low redshift contain a 
significant amount of baryons. 

A second Nevm detection (log[o(NNcviii/cm"^) = 13.98, = 
47 mA) at z = 0.326 was reported by Narayanan et al. (2009). This 
absorber was identified in a FUSE spectrum of the QSO 3C 263 
(Zem = 0.646) with S/N = 5-10. This group found that the ionisation 
state of the gas phase giving rise to the observed Ne vni absorption 
could not be explained by photo-ionisation. Rather, this absorber 
appears to be multi-phase, with the Ne vm absorption arising in 
gas at r ~ 10* K. These authors provided the first estimate of 
the baryonic content of Nevm systems, fi;(,(Neviii) = 0.00263, 
which corresponds to roughly 6 per cent of the total baryons in 
the Universe (assuming = 0.0418). This absorber was later 
confirmed in a COS spectrum by Narayanan et al. (2012). They re- 
analysed the physical conditions in this O vi- Ne vm absorbing gas, 
and found that is is best modelled with gas in CIE at T = 5.2 x 10^ K 
with [Ne/H] = -0.12+^;}2 (adopting (nNe/«H)o = -4.07). These 
authors noted for the first time the surprisingly small mmiber of 
Ne vm detections at low redshift'. This fact led to the interpretation 
that the Ne vm absorbers have a physical origin different from the 
filamentary structures of the highly ionised gas outside the virial 

' For comparison, there are more than 100 O vi detections at low redshift 
reported in the literature to date (see e.g. Tilton et al. 2012). 



boundaries of galaxies. This group suggested that the observed 
Ne vm absorption originates in the extended regions around galaxies, 
i.e. the circum-galactic medium (COM), which is consistent with 
the relatively high metal abundances of the Ne vm absorbers. 

Narayanan et al. (201 1) next reported on the detection of ab- 
sorption by Ne vm (and associated O vi) in an intervening, multi- 
phase system at z = 0.495 in a COS spectrum (S/N=15 - 20) along 
the sightline to the source PKS 0405-123. Using the apparent optical 
depth method (AOD; Savage & Sembach 1991), they measured 
logiQ(NNeviii/cm"^) = 13.96 and = 45 mA. Allowing for the ob- 
served O VI to have a significant contribution from the photo-ionised, 
cooler gas phase, the ionisation state of the warmer, O vi- Ne vm ab- 
sorbing phase can be explained with a model that takes into account 
both photo- and collisional ionisation. Adopting a fixed temperature 
T = 4.7 X 10^ K and [Ne/H] = -0.3 (relative to (;jNe/nH)o = -4.07), 
Narayanan et al. (201 1) estimated a gas density Wh ~ 8 X lO^'^cm"^, 
and a linear size of the absorber of roughly 20 kpc. 

Another prominent example is the system in a COS spectrum 
of the QSO PG 1206+459 reported by Tripp et al. (201 1). This 
group detected a strong, complex Ne vm absorber (and absorption 
by a series of other high and low ionisation states of other ele- 
ments: Hi; Mgn, Mgx; Nii, Niii, Niv, Nv; Om, Oiv; Sm, Srv, 
S v) at z = 0.927, consisting of nine individual components with a 
total column density logio(NNcviii/cm"^) « 14.9. They found strong 
evidence that the Ne vm absorption is produced by the highly en- 
riched, kinematically complex plasma in a strong outflow from a 
post-starburst galaxy (L = 1.8 U) at z = 0.927 and a projected 
impact parameter p = 68 kpc from the sightline. They estimated 
that the warm-hot (r > 10^ K) gas phase contains 10 - 150 times 
more mass than the cool (J < 10^ K) photo-ionised gas phase, with 
individual absorbing components harbouring 10* Mg - 4 x 10'" Mg. 

Finally, three intervening systems at z = 0.684, 0.705, and 
0.725 along a single sightline to the QSO PG 1 148+549 containing 
Nevm (and Ovi) were reported by Meiring et al. (2012), These 
absorbers were detected in a COS spectrum with S/N = 20 - 
40. Two of the absorbers (z = 0.684, z = 0.705) appear to be 
two-component systems, with total log|(,(N>jevra/cm"-) = 13.95 
and W,. a; 43 mA, and total log|(,(NNeviii/cm"^) = 13.86 and 
Wr ~ 35 mA, respectively. The system at z = 0.725 shows a very 
simple kinematic structure (probably single-component) with 
logio(NNevin/cm~^) = 13.81 and Wr ~ 33 mA. All three Ovi - 
Ne vm absorbers can be explained as being produced in gas in CIE 
with a temperature around log|(,(r/ K) = 5.7, and a metallicity (i.e. 
oxygen abundance) [O/H] > -0.5. Assuming the Ovi - Nevm 
absorption is produced in the CGM of galaxies, and using the 
observed absorber redshift density and the galaxy luminosity 
function, Meiring et al. inferred an average size of 70 - 150 kpc 
for typical Ovi - Nevm absorbers at z = 0.7. This group further 
estimated that the density of gas as traced by these absorbers is 
Hj, = 0.002, or roughly 4 per cent of the baryon density (assuming 
Qb = 0.0418). 



APPENDIX B: SYNTHETIC SPECTRA 

The Nevm line statistics (Section 2.1.2) and the Nevm- and Hi- 
optical depth weighted quantities used throughout this work have 
been extracted from a set of spectra along 5000 random sightlines 
drawn from our simulation at z = 0.5, filling a total redshift path 
Az = 212, corresponding to an absorption path length Ax = 387. 
For each sightHne we have computed a pair of synthetic spectra. 
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Figure Al. Left: Distribution of Doppler parameters of the Ne viii components. The dashed vertical line indicate in each case the resolution limit, i.e. the minimum 
adopted fc-value. The blue line show the distribution of components identified in spectra that have been convolved with a Gaussian LSF {FWHM =17 km s" ' ), 
while the red line shows the distribution of components identified in spectra that have not been convolved. Centre: Comparison of the total line width, Z^ncviii, to 
the thermal line width, bj (equation 9) of the single-component Ne viii components identified in high-S/N spectra in our simulation at z = 0.5. Clearly, the 
intrinsic width of the absorption lines is mainly set by temperature of the absorbing gas. However, convolution with an instrumental LSF results in absorption 
lines that are broader with respect to the pure thermal line width. Note that roughly half (45 per cent) and one third (35 per cent) of the Ne vni absorbers appear 
as single-component features in the convolved and non-convolved spectra, respectively. Right: Rest-frame equivalent width and column density for our sample 
of Ne VIII absorbing components identified in 5000 spectra along random sightline drawn from our simulation at z = 0.5; each of the black curves shows a 
theoretical curve-of-growth (CoG) adopting foNcviii = 50 km s"' , 18 km s"' , and 4 km s"' (from top to bottom). Meaning of colours identical to the left panel. 
Note that there is a linear relationship between equivalent width and column density for all Ne viii components. Since equivalent width is conserved, convolution 
with an instrumental LSF does not affect the measured column density. 



with each spectrum containing absorption by Hi Lya or Neviii 
770.41 A only. These spectra have been generated using the package 
SPECWIZARD written by Joop Schaye, Craig M. Booth, and Tom 
Theuns, following the method described in detail in Tepper-Garcia 
et al. (201 1, their Section 3.1). We have convolved our spectra with 
a Gaussian line-spread function (LSF) with a FWHM =17 km s"' , 
and have re-sampled our spectra onto pixels of 3kms"' in size, 
which toughly match the instrumental properties of COS. Finally, 
we have added Gaussian noise with a flux dependent root-mean- 
square (rms) amplitude given by (S/N)"' F{v), where F{v) is the flux 
as a function of velocity along the sightline, and S/N is the adopted 
signal-to-noise ratio; we adopt S/N=1000 to allow for the detection 
of most of the Nevin in our simulation (see below and Section 
2.2). Throughout this work, we shall refer to these as 'high-S/N, 
COS-like' spectra. 

We fit our synthetic spectra using a significantly modified ver- 
sion of AUTOVP (Dave et al. 1997) following the method detailed in 
Tepper-Garcia et al. (2012, their Appendix A). Briefly, Ne viii ab- 
sorption lines are identified in regions with an (integrated) rest-frame 
equivalent width W,- > N crv/^. , where criy,. is the uncertainty in the 
equivalent width, integrated over n pixels, and N is the significance 
level. We adopt n = 25 (equivalent to 75kms"' or approximately 
four times our adopted spectral resolution) and'" N = 5. With these 
values, the significance value chosen translates into a formal, rest- 
frame equivalent width limit W,- ~ 193 (S/N)"' mA. This implies 
a formal completeness limit W,- a 0.2 mA or log|Q(NNeviii/cm"^) » 
11.5 for our spectra with S/N=1000. To ensure the detection 
of all relevant features, we impose a minimum column density 
logjo(NNeviii/cm"') = 11.0, which is a factor 3 below the formal 
completeness limit. We note that the effective completeness limit is 
somewhat higher, log(()(NNeviii/cm"-) a 11.7 or W,- ~ 0.3 mA. 

Each component is fitted assuming a Gaussian profile," and the 
parameters (velocity centroid vq, column density Nncviii, and Dop- 



pler parameter ^incvoi) of all lines identified in each region are sim- 
ultaneously and iteratively adjusted using the Levenberg-Marquardt 
algorithm (Levenberg 1944; Marquardt 1963) as implemented in 
Press et al. (1992) until the reduced ;^'^-value, i.e. the;^'^-value di- 
vided by the degrees of freedom, is below = 1.25. During each 
iteration, lines with relative errors in N^e vni or VIII larger than 50 
per cent are discarded. Note that each region is fitted with the least 
number of absorption components consistent with the condition 
< X'^f^- Fitting of our 5000 simulated spectra yields a total of 8665 
Ne viii components, with roughly half (45 per cent) of the identified 
Ne viii absorbers being single-component features. The total cumu- 
lative line frequency per unit redshift is dN/dz(Wr > 0.3 mA) a 41. 

To investigate the effect of the instrumental broadening on the 
line properties, we fit a second set of spectra which are identical 
to our COS-like spectra except that these new spectra are not con- 
volved with an instrumental LSF. In Figure Al we compare the 
line parameters of the two Ne viii absorption lines samples obtained 
from these two different sets. The left panel demonstrates that con- 
volution with a Gaussian LSF artificially increases the line width 
by a significant amount. Therefore, while the intrinsic line width is 
a on average good indicator of the gas temperature (central panel, 
inset), the width of lines identified in convolved spectra can at most 
provide an upper limit to the gas temperature. The right panel shows 
the relationship between rest-frame equivalent width and column 
density (i.e. the curve-of-growth). Clearly, the equivalent width of 
the majority of the Ne viii absorbers depends linearly on the Ne viii 
column density. Since the rest-frame equivalent width is conserved 
under a convolution, the column density of the lines is conserved as 
well. 



The Ne viii detections reported to date in the literature have a significance 
No- > 4a-. 

" The column density of the lines are low enough that the difference 
between a Gaussian and a Voigt-Profile is negligible. 
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